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I. INTRODUCTION 

Vortices have been a source of fascination since the 
works of Empedocles, Aristotle and Descartes, who tried 
to explain the formation of the Earth, its gravity and the 
dynamics of the solar system as due to primordial cos- 
mic vortices. Many interesting problems related to vor- 
tices are open in different fields such as fluid mechanics, 
superconductivity, superfluidity, light propagation, Bose- 
Einstein condensation (BEC), cosmology, biosciences, or 
solid state physics [J, i, i, i, (i, . 

In wave mechanics a vortex is a screw phase disloca- 
tion, or defect J7|], where the amplitude of the field van- 
ishes. The phase around the singularity has an integer 
number of windings, ^, which plays the role of an an- 
gular momentum. For symmetric systems, this number 
is a conserved quantity and governs the interactions be- 
tween vortices as if they were endowed with electrostatic 
charges. Thus, £ is usually called the topological charge 
of the defect. 

Angular momentum is conserved in a quantum sys- 
tem with 0(2) rotational symmetry. If we consider a 
state with well-defined angular momentum £ G Z, i.e., 
an eigenfunction of the angular momentum operator at a 
given time to, its evolution will preserve the value of £. In 
a system possessing a discrete point-symmetry (described 
by the Cn and Cnv groups) the angular momentum is no 
longer conserved. However, in this case one can define 
another quantity m G Z, the Bloch or angular pseudo- 
momentum, which is conserved under time evolution [8] . 
The angular pseudo-momentum m plays then the role of 
^ in a system with discrete rotational symmetry. From 
the group theoretical point of view, the angular momenta 
and pseudo- moment a, i and m, are also the indices of the 
2D irreducible representations of 0(2) and Cn, respec- 
tively [9,|lO,[ii|. Unlike ^, the values of m are limited by 
the order of the point-symmetry group Cn- 

The existence of vortices is one of the signatures 
of superfluidity and this is why in the field of Bose- 
Einstein condensation, they have attracted so much in- 
terest. They can be generated in rotating traps [Bl, @ 



or by phase imprinting methods [l2l, K§\ . The later pro- 
cedure allows to generate only multiply charged vortices 
with topological charges m = 2, 4. 

In the last years there has been an enormous interest 
on the applications of group theory to study the proper- 
ties of defects in media with discrete symmetries includ- 
ing photonic crystals, periodic potentials with discrete 
symmetries, etc jl, ll, |lO|, 0, [l6| . 

In this paper we explore the application of group the- 
ory to control the topological charge of vortices in Bose- 
Einstein condensates by using external potentials with 
discrete rotational symmetry. To do so we propose a sim- 
ple setup based on a non-periodic potential with discrete 
rotational symmetry which will allow us to perform many 
operations with the vortex charges depending on the ini- 
tial charge and the potential symmetry order. While the 
vortex transmutation has been previously explored in po- 
tentials with broken symmetries [18] and in the context 
of photonic lattices [14], other operations to be proposed 
here have not been studied before. Our proposal is sim- 
pler to implement than the 2D lattice type potentials 
proposed in the framework of photonic lattices [15] and 
is easier to reconfigure. 

We will also show how starting from multiply charged 
vortices such as the ones which can be_generated in atom 
chips by phase-imprinting methods [UllH one can gen- 
erate different types of vortices by choosing an appropri- 
ate control potential. 

Our plan is as follows: In Sec. [Ill we present the theory 
in which our methodology is based. In Sec. [Ill we present 
several examples. First we present our specific setup in 
Sec. nil Al and then discuss several phenomena which 
can be achieved such as: topological charge erasing (Sec. 
nil Bp . and two examples of topological charge inversion 
from m = 2tom = —1 (Sec. IIII C\i and m = 4 to 
m = -l rSec. [HTDl) . 
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II. THEORY 

In this paper we will consider a Bose-Einstein conden- 
sate with tight confinement along an specific direction (z) 
leading to a quasi-two dimensional Bose-Einstein conden- 
sate. Under the effect of an additional external potential 
y(x) this system is ruled in the mean field limit by an 
effective Gross-Pitaevskii equation [17^ lIS] 

idMx) = -^Ax^(x) + V{^)^{x) + 5|^(x)| V(x), (1) 

in dimensionless units and where ^ is a measure of the 
effective nonlinearity. 

The symmetry of the potential induces strict restric- 
tions on the vorticity of stationary solutions of Eq.([T]). 
In order to illustrate this statement, let us consider the 
equation for stationary solutions, which is a nonlinear 
eigenvalue equation of the following type: 

if(V'(x),x)V(x)=Mx), (2) 

where the nonlinear Hamiltonian operator depends on 
the field itself, i.e., (^(x), x) = -|Ax+V(x)+^/|V^(x)p 
The nonlinear solution ip can be considered as a self- 
consistent solution. It appears as an eigenstate of the 
Hamiltonian but, at the same time, it defines the Hamil- 
tonian operator itself since H depends on it through the 
nonlinear term ^IV^p. 

One can establish necessary conditions for the exis- 
tence of stationary symmetric solutions of Eq.([2]) based 
on the symmetry properties of the potential V{'x.). They 
will exhibit special features due to the nonlinear nature 
of the Hamiltonian operator H{2p). We can summarize 
the properties of these solutions in a single statement: 
(i) if the potential V is invariant under a point symme- 
try group that we refer to generically as Q — i.e., Q de- 
scribes finite two-dimensional 2i\ jN rotations around an 
axis (Cat group) and specular refiections (Cnv group) — 
and (ii) if we search for symmetric solutions fulfilling the 
condition |?/^(Gx)p = |?/^(x)p (where G is any element of 
the group Q\ then the solution il) must belong to some 
representation D^{Q) of the symmetry group Q or to 
some of their subgroups C Q. Since this statement is 
rather mathematical, it is convenient to analyze it in the 
light of the symmetry properties of the nonlinear Hamil- 
tonian H(ip). First of all, let us recall that if is a given 
stationary solution = V^soi satisfying Eq.([2]) then V^soi 
plays two different roles in this nonlinear eigenvalue equa- 
tion. On the one hand, ?/^soi defines the Hamiltonian op- 
erator i^(V^soi) whereas, on the other hand, it appears as 
an eigenfunction of the same operator. This two- fold role 
has profound implications on the allowed functional form 
of the solution. The first consequence of the explicit, and 
specific, dependence of the nonlinear Hamiltonian on ip^^i 
is that H(ipQo\) inherits the symmetry of the potential. 
This is a consequence of assumptions (i) and (ii). Since 
the Laplacian operator is invariant under any type of ro- 
tation Agx = Ax and, on the other hand, V(Gx) = V{^) 



and |?/;soi(Gx)p = l^/^soKx)!^ because of the previous as- 
sumptions, the nonlinear Hamiltonian evaluated at ?/^soi 
is automatically invariant under the symmetry group 
g: H{^{G^),G^) = -|AGx + nGx) +^|^(Gx)|2 = 
-^Ax + l^(x)+^|V^(x)p = (V^(x),x). The fact that the 
Hamiltonian operator evaluated at ?/^soi is invariant un- 
der the group Q implies that H{ipso\) commutes with all 
elements of this group, i.e., [i^(?/^soi), G] = VG G G and, 
therefore, according to standard quantum mechanics ar- 
guments, all its eigenfunctions must belong to the dif- 
ferent representations of the symmetry group Q. At this 
point, we have considered i/jsoi only in its role of generator 
of the nonlinear Hamiltonian operator. However, in its 
second role ?/^soi must also appear as an eigenfunction of 
H{iIjso\)' Therefore, i/jsoi must belong to some of the rep- 
resentations of i^(V^soi), i-e., V^soi ^ D^{Q), where in- 
dicates the representation characterized by the represen- 
tation index m and the degeneracy index u. The degener- 
acy of the representations of the point symmetry groups 
Gat and Cnv is either one or two [lo|. One dimensional 
representations automatically fulfill the symmetry condi- 
tion for the amplitude |?/^(Gx)p = |?/^(x)p assumed in (ii) 
since they do not transform, except for a sign, under the 
action of a finite rotation: GV^(x) = tjj{Gx.) = ±?/^(x). 
The requirement of amplitude invar iance is, however, 
trickier for two-dimensional representations. Let us see 
next why. A natural basis for a two-dimensional repre- 
sentation of a point symmetry group is that formed by 
the eigenvectors of the group operator in such represen- 
tation, i.e., that given by the two functions (i/jm^i^m) 
fulfilling G'?/^rn(x) = ilJm{Gx) = e^?/;^(x) and its com- 
plex conjugate, m being the representation index and 
the eigenvalue being given by e = e*^^/^. Since a group 
rotation operator acting on this representation will be 
represented by a 2 x 2 matrix this two vectors provides 
the basis where this matrix is diagonal. However, the 
more general form of a function belonging to the m- 
representation will be given by a linear combination of 
ipm and ip^. This general solution has a problem with 
respect to the requirement of invariance for the ampli- 
tude assumed in (ii). Indeed, in the most general case in 
which we consider arbitrary coefficients, this linear com- 
bination does not verify the previous condition and thus 
such a function would be excluded as a solution of the 
problem enjoying symmetry under the full group Q (how- 
ever, it can fulfill the condition for a subgroup Q' C 
see Ref.[9]). In this two-dimensional subspace the only 
two functions fulfilling the condition are the eigenvectors 
of the G operator i/j^ and ip^. This is so because the 
eigenvalue is a pure phase number and therefore with 
unit modulus. Consequently, only tjjm and tjj^ can be 
nonlinear solutions of Eq.([2]). This fact is in remarkable 
contrast with respect to the linear case in which all linear 
combinations belonging to the representation appear as 
solutions of the eigenvalue equation. Notice that the tjjm 
and t/j^ complex solutions represent a vortex- ant ivortex 
soliton pair of topological charge m and —m respectively. 
In order to interpret this result on a more physical basis it 
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is convenient to re-write the action of a discrete rotation 
of order N on these functions using polar coordinates: 

G4>^{r, 6) = 4>^{r, 6 + = e^2™/^^^(r, 9) (3) 

and its complex conjugate. As it was recognized in 
Ref. ([!]), this condition is identical to that fulfilled by 
one-dimensional Bloch modes but by substituting a stan- 
dard non-compact spatial coordinate by the compact an- 
gular one 6. For this reason the most general form for 
such solutions are that corresponding to angular Bloch 
modes: 

^m{r,e)=e'^'fm{r,e) (4) 

where fm is an invariant function under 2it/N rotations 
and m plays the role of "angular" pseudo-momentum sat- 
isfying the following constraints: \m\ < N/2 (if N is 
even) and |m| < (TV — l)/2 (if is odd). Since m is 
also a representation index, these constraints can be ex- 
plained using group theory arguments [11]. However, it 
is physically more intuitive to interpret them by means 
of the properties of the angular pseudo-momentum in an 
equivalent Bloch problem in terms of the angular vari- 
able 0. The symmetry condition of the potential under 
A/'th-order rotations — V{r^ + 27r/N) = V(r, 0) — is un- 
derstood as a periodicity condition in the angular vari- 
able 6 with period given by a = 2t:/N. As we have seen, 
this periodicity property is inherited by the full Hamil- 
tonian i^(V^) when fulfills the symmetry condition 
|V^(r,6> + 27r/A')P = |V^(r,6>)p (assumption (ii)). Since the 
full Hamiltonian is periodic it is then natural that their 
solutions transform as in Eq.Q and they present the an- 
gular Bloch form (|4|). The role of m is therefore that cor- 
responding to the conjugate variable of the periodic vari- 
able 6>, hence its name of "angular" pseudo-momentum. 
The angular nature of forces m to be an integer in order 
to preserve the single- valuedness of the wavefunction ?/^. 
Energy is a periodic function of pseudo-momentum with 
periodicity given by the so-called reciprocal lattice vector 
27r/a = A", which determines the existence of equivalent 
zones in pseudo-momentum space called Brillouin zones. 
In this framework the constraints in the values of m arise 
from the fact that all possible solutions can be reduced 
to those existing in the first Brillouin zone in pseudo- 
momentum space, i.e., those fulfilling \m\ <T:/a = N/2, 
which is equivalent to the aforementioned constraints. As 
an example, in FigH] it is represented the angular Bloch 
structure of energy eigenstates of two systems owning 
discrete symmetry of 4th and 5th order. Notice that the 
period of the Brillouin zone — i.e., the reciprocal lattice 
vector — is equal to the order of symmetry (A" = 4, 5) 
and that the allowed angular pseudo-momenta fullfill the 
condition m < 2 corresponding to the first Brillouin 
zone. In this scenario it is rather intuitive to under- 
stand what happens to a solution carrying angular mo- 
mentum / evolving in a medium characterized by full ro- 
tational invar iance — given by the 0(2) group — when 
we suddenly switch on a potential that breaks full rota- 
tional symmetry into a discrete-symmetry of A^th order 



— given by the Cn group. In terms of the angular vari- 
able the solution with angular momentum / behaves 
like a sort of "angular" plane wave e^^^ fi{r) since the 
amplitude fi(r) is angle-independent. The angular mo- 
mentum / plays the role of ordinary (discretized) mo- 
mentum. The problem is thus equivalent to that of a 
plane wave propagating in a constant potential — in 0, 
since U{r) = V(r) + gfi{r) is angle independent — that, 
at some specific moment to, feels the presence of a peri- 
odic potential U{r,0 -\- a) = U (r, 0). At to the wavefunc- 
tion corresponding to the initial angular plane wave must 
excite the spectrum of the nonlineal operator H{ilj{to)). 
However, the fuh potential U{r, 0; to) = V(r, 0) + gfiir) 
is no longer 0(2)-invariant but CAr-invariant due to the 
breaking of continuous symmetry by the appearance of 
the discrete-symmetry potential V{r, 0) at to and whose 
eigenstates are angular Bloch modes of the form (|4]). 
This fact implies that there must be a matching between 
the input angular momentum and the angular pseudo- 
momenta of allowed angular Bloch states after the Cn 
interaction is switched on. In other words, the initial 
angular momentum / carried by the angular plane wave 
e^^^ fi{r) must match the angular pseudo- momentum m 
of some angular Bloch state. If |/| < A'/2 the initial an- 
gular momentum can always match an angular pseudo- 
momentum in the first Brillouin zone and, consequently, 
/ = m. However, if |/| > A^/2 we excite angular Bloch 
states in higher-order Brillouin zones which, on the other 
hand, we know are equivalent to those in the first Bril- 
louin zone. This equivalence is given by periodicity in 
pseudo-momentum space which is provided by the recip- 
rocal lattice vector 27r/a = A^ in such a way two pseudo- 
momenta are equivalent if they differ from each other 
by a multiple of this vector. This establishes the de- 
sired matching condition for angular pseudo-momentum 
in terms of the initial angular momentum /: 

l-m = kN, keZ. (5) 

This matching condition can be easily visualized in 
FigUl for the A" = 4 and A" = 5 cases. When A" = 4 
and / < 2 we can always find Bloch states in the first 
Brillouin zone with the same value of the angular pseudo- 
momentum m = I. This is the case of the input angular 
momentum / = 2 represented in Figdfa). For N = 5 the 
situation is similar, i.e., we excite angular Bloch states 
with m = / as long as / < 2 (see Figdjb)). However, 
when this condition is not fulfilled (/ > 2) we excite 
an equivalent Bloch state with different angular pseudo- 
momentum in the first Brillouin zone. Let us choose, for 
example, the input angular momentum to be / = 3 in the 
N = 5 case. We immediately see that we are exciting 
an equivalent Bloch state with different angular pseudo- 
momentum in the first Brillouin zone given by m = —2, 
in agreement with the matching condition ([5]). 

The previous considerations refer to angular pseudo- 
momentum and thus they are related to the way the 
wavefunction transform under discrete rotations; i.e., in 
the way described by Eq.(j3|). Using analogous arguments 
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FIG. 1: [Color online] Angular Bloch structure of solutions 
for a system with rotational symmetry of order: (a) A/" = 4 
and (b) N = 5. Red and green circles symbolize solutions 
that belong either to one dimensional or two dimensional ir- 
reducible representations. Red arrows represent input angu- 
lar momentum while blue arrows represent output angular 
pseudo- momentum. Some examples of matching conditions 
for angular pseudo- momentum are indicated. 



directly from inspection of the bloch functional form (|4]). 
In fact, if N is even, the dimension of the representation 
is one for m = or | m \= while it is of dimension 
two for I m 1= 1, . . . , ^ — 1. And if N is odd, the di- 
mension of representation is one for m = 0, while it is 
two for I m 1= 1, . . . , Since one dimensional repre- 

sentations do not transform under rotations except for a 
sign, they are necessarily real. On the other hand, two 
dimensional representations present a non trivial phase 
structure. Thereby, only solutions belonging to two di- 
mensional irreducible representations can be considered 
as vortex solutions. Notice that, as stated previously, the 
base for these two dimensional irreducible representations 
is the complex pair ?/^^ and t/j^ of vortex- ant ivortex soli- 
ton solutions. 

According to these statments and to the matching con- 
dition ([5]) there are different kinds of vorticity transfor- 
mations between 0(2) and Cnv media. For / in the first 
Brillouin zone no transformation is produced. For / out- 
side the first Brillouin zone there are two kinds of trans- 
formation: i) if the value of m that corresponds to the 
input / according to the matching condition ([5]) is such 
that the solution belongs to a irreducible representation 
of dimension one, the phenomenon is called charge eras- 
ing and ii) if m is such that the related representation is 
a two dimensional one, the phenomenon is called vortex 
transmutation. Moreover, there are two kinds of vortex 
transmutation phenomenon: i) charge downconversion if 
m has the same sign that / and ii) charge inversion if m 
has different sign that I. For example, in in FigUJa) a 
case of charge erasing is presented. On the other hand, 
the Figlljb) represents a charge inversion. 

This remarkable fact can be used to exploit the angular 
pseudo-momentum matching condition ([5]) as a way to 
control the topological charge of matter wave vortices, as 
it will demonstrated extensively in the next sections. 



III. APPLICATIONS 



as those appearing in Ref.[8], it can be proven that the 
angular pseudo-momentum is conserved during time evo- 
lution. This means that if we start at the initial time to 
with a solution fulfilling Eq. (|3|) , this solution will preserve 
this property during evolution provided no new alteration 
in the symmetry of the potential is given. In other words, 
if the initial angular pseudo-momentum is m at to it will 
remain the same at any other time t > to- This property 
has been numerically checked in optical vortex transmu- 
tation phenomena [ij]. When, besides having a solution 
with angular pseudo- momentum m, we are dealing with 
a function characterized by a single-phase singularity the 
angular pseudo-momentum m becomes identical to the 
topological charge [11, 14]. 

It is remarkable that the irreducible representations in 
Cnv are either one or two dimensional. This statement 
is based in group theory arguments or can be deduced 



A. Practical configuration 

In order to work with the minimal configuration show- 
ing the phenomena to be described here we have made 
the simplifications of taking the linear limit ^ = and 
chosing a simple potential V{x^y) obtained as the super- 
position of N gaussian functions of the form 

N-l 

i=o 

(6) 

with {xj^yj) = d (cos(27rj/A^), sin(27rj/A/')) This type of 
potential can be obtained physically by using a set of 
laser beams generating a set of optical dipole traps for 
the Bose-Einstein condensate. 

We have checked that our results are essentially inde- 
pendent of the nonlinearity and that they remain valid 
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FIG. 2: [Color online] Erasing of a vortex of topological charge 
m = 2 due to the effect of the discrete symmetry of the po- 
tential ([6]) with N = A^d — 4,Vb = —0.6. Shown are pseudo- 
color plots of the amplitude |'0(x,y,t)|^ (a,b,c,d) and phase 
arg('0(x, t)) (e,f) for different times indicated on the sub- 
plots. The locations of the potential minima according to Eq. 
(|6]) are indicated with small white circles. 



B. Charge erasing 

As a first example of vorticity control by discrete sym- 
metries we consider the evolution of an initial configura- 
tion of the form 

^{x, y) = {x^ iyf exp [{-x^ - y^)/^] (7) 

i.e. a vortex of topological charge 1 = 2. This initial 
wavefunction will be subject to a potential with fourth 
order symmetry corresponding to Eq. (|6]) with N = 4. 
The other potential parameters are taken to be Vb = —0.6 
and d = A. 

The evolution of this initial datum is shown in Fig. [2l 
Where some typical features are seen. First of all the 
amplitude of the wavefunction [Fig. [2ja,b,c,d)] becomes 
localized on the four potential wells after a transient in 
which some radiation is emitted. To achieve this effect of 
localization it is convenient to adjust the potential min- 
ima to be close to the initial maximum density so that 
the fraction of the initial number of particles which is 
retained into the system is maximized. 

At the same time the phase [Fig. [2fe,f)] experiences a 
complicated evolution from the initial configuration with 
/ = 2, to a final configuration with two nodal lines with 
X = and ^ = (which are the symmetry axes of this 
potential) and constant phases across the four regions in 
which the space is separated by them, which correspond 
to a solution with m = 2 that belongs to a representation 
of dimension one. Let us obtain the theoretical prediction 
for this case. We must take /c = in expression ([5]), 
since the input vortex is in the first Brillouin zone, and 
consequently, m = 2 — 0-4 = 2, in agreement with the 
numerical result. 



for the more complicated case of a periodic potential of 
the required symmetry. However, the choice of the po- 
tential as given by Eq. ([6]) is not only simpler but allows 
for more freedom in the selection of the symmetry since 
we can take e.g. = 5 to obtain a discrete symmetry of 
fifth order, something which is not possible with lattice 
potentials. 

It is also very simple to change between configurations 
with different symmetry orders (A") even dynamically by 
just adding or eliminating laser beams. 

In what follows we will solve Eq. ([1]) with V(x, ?/) 
given by Eq. (|6|) and initial data of (multicharged) vor- 
tex type. All simulations to be presented in this paper 
have been done using a split-step method where the spa- 
tial derivatives are computed by using a pseudoespectral 
formula based on trigonometric polynomials. In all cases 
to be presented here we have chosen At = 0.025 and 
the simulation region [—20, 20] x [—20, 20]. In the figures 
the spatial region shown is [—10,10] x [—10,10]. The 
outgoing radiation is elliminated by absorbing boundary 
conditions such as the ones implemented in Ref. [19j . 



C. Generation of singly charged vortices from 
doubly charged ones: charge inversion 

An interesting problem is how to process the doubly 
charged vortices obtained in atom chips by phase im- 
printing methods [l2|, [l^ to obtain vortices with charge 
m = —1. 

As in the previous examples, discrete symmetries can 
be helpful to accomplish this task. Let us take now a 
potential with A" = 3 and a doubly charged vortex as 
initial datum as in Eq. ([7j). Let us choose the potential 
parameters as d = 3.2, and Vb = —0.6. 

The evolution of this configuration is shown in Fig. [3] 
where the inversion of the central vortex from / = 2 to 
m = — 1 is clearly seen in the phase plots [Fig. [3l^b,d,f,h)]. 
In this case, to obtain the theoretical prediction, we must 
take /c = 1, since / = 2 falls in the second Brillouin zone 
for A" = 3. Consequently, according to the matching con- 
dition ([5]),m = 2 — 1-3=- 1. It can be seen how initially 
a central vortex with charge m = — 1 coexists with other 
vortices located around it. However, as time becomes 
larger these vortices slowly spiral out of the system and 
are not visible in the region of interest for t ^ 80. 
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FIG. 4: [Color online] Vorticity inversion from m = 4 to 
m = — 1 due to the effect of the discrete symmetry of the 
potential (|6|) with N — 5,d — 3.2, Vo — —2. Shown are pseu- 
docolor plots of the amplitude 2/, t)|^ (a,c,e) and phase 
arg('0(x, ^, t)) (b,d,f) for different times indicated on the sub- 
plots. The locations of the potential minima according to Eq. 
(|6|) are indicated with small white circles. 



FIG. 3: [Color online] Vorticity inversion due to the effect of 
the discrete symmetry of the potential JG]) with N — — 
3.2, Vb = —0.6. Shown are pseudocolor plots of the ampli- 
tude 2/, t)|^ (a,c,e,g) and phase arg('0(x, 2/, t)) (b,d,f,h) 
for different times indicated on the subplots. The locations of 
the potential minima according to Eq. (|6]) are indicated with 
small white circles. 



D. Obtention of a vortex of charge m — —1 from a 
vortex with charge / = 4 in a system with symmetry 
= 5: charge inversion 

As a final example we show how to obtain a vortex with 
charge m = — 1 from a vortex with topological charge 
/ = 4. To do so we must use a system with symmetry 

= 5. Our initial wavefunction is given by 



The transmutation phenomenon has been studied pre- 
viously in Refs. [flllll- However, in Ref. [18] the asym- 
metric trap used to induce the phenomenon leads to a 
recurrence in which the topological charge periodically 
oscillates between +1 and -1, while in this setup we get 
m = — 1 for any time larger than the transient time in 
which the system becomes stabilized. With respect to the 
proposal of Ref. [14] the scheme presented here has the 
advantage of introducing less radiation and being cleaner 
because of the finite range of the potential. 



^{x, y) = {x^ iyf exp [(-x^ - y')l2\ , (8) 

and the parameters in the potential are taken as = —2 
and d — 3.2. 

As in the previous examples we observe emission of ra- 
tiation (in this case stronger due to the increasing energy 
of the initial condition which leads to a faster scape from 
the origin) and after some time, the desired structure is 
observed. Again, /c = 1 since / = 4 falls in the second 
Brillouin zone for A/" = 5. Thereby, according to ([5]), 
m = 4—1-5 = —1, in agreement with the numerical 
solution. 
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IV. CONCLUSIONS AND DISCUSSION 

In this paper we have studied the dynamics of multiply- 
charged vortices under the action of potentials with dis- 
crete symmetries. We have shown how the symmetry or- 
der of the potential can be used to manipulate the topo- 
logical charge of the vortex in order to obtain a desired 
(lower) value of the vorticity. Specifically we have pro- 
posed a setup of gaussian traps which can be configured 
to have the appropriate symmetry and more fiexibility 
than the lattice type potentials previously used. 

As applications we have studied: topological charge 
erasing, charge inversion from / = 2 to m = — 1 and 
topological charge conversion from / = 4 to m = — 1 
The last two examples are applicable to the controlled 
generation of singly charged vortices from the output of 
the phase imprinting method using matter-wave chips 
developed in Ref. |12l |. 



The vortex charge control and erasing properties which 
can be achieved in this system open new possibilites for 
the control of quantum matter. Moreover, although we 
have choosen an specific model given by Eqs. ([T]) and 
(j6|) our arguments developed in Sec. [TTl are completely 
general and depend only on the discrete symmetry prop- 
erties of the system. Thus, they can be extended beyond 
the specific form of the potential choosen and even the 
Gross-Pitaevskii equations used in this paper. 

Acknowledgments 

This work has been partially supported by grants 
BFM2003-02832, FIS2004-20188-E, FIS2005-01189, and 
FIS2006-04190 (Ministerio de Educacion y Ciencia, Spa- 
in), PAI05-001 (Consejeria de Educacion y Ciencia de la 
Junta de Comunidades de Castilla-La Mancha). 



[1] H. J. Lugt, Vortex Flow in Nature and Technology, 

Krieger, Malabar, FL, (1995). 
[2] L. M. Pismen, Vortices in Nonlinear Fields, Clarendon, 

Oxford, UK, (1999). 
[3] F. Sols, Physica C 369, 125 (2001). 

[4] M. R. Matthews, B. P. Anderson, P. C. Haljan, D. S. 
Hall, C. E. Wieman, E. A. Cornell, Phys. Rev. Lett. 83, 
2498 (1999). 

[5] K. W. Madison, F. Chevy, W. Wohlleben, and J. Dal- 

ibard, Phys. Rev. Lett. 84, 806 (2000); F. Chevy, K. W. 

Madison, J. Dalibard, Phys. Rev. Lett. 85, 2223 (2000). 
[6] J.R. Abo-Shaer, C. Raman, J.M. Vogels, W. Ketterle, 

Science 292, 476 (2001). 
[7] J. F. Nye, M. V. Berry, Proc. R. Soc. London A 336, 

165-190 (1974). 
[8] A. Ferrando, Phys. Rev. E 72, 036612 (2005). 
[9] A. Ferrando, M. Zacares, P. Fernandez de Cordoba, J. A. 

Monsoriu, and P. Andres, Opt. Express 13, 1072 (2005). 
[10] M. Hamermesh, Group Theory and its Application to 

Physical Problems, Addison- Wesley Series in Physics 

(Addison- Wesley, Reading, MA, 1964), 1st ed. 
[11] A. Ferrando, M. Zacares, and M. A. Garcia-March, Phys. 



Rev. Lett. 95, 043901 (2005). 
[12] A. E. Leanhardt , A. Gorlitz, A. P. Chikkatur, D. Kielpin- 

ski, Y. Shin, D. E. Pritchard, and W. Ketterle, Phys. 

Rev. Lett. 89, 190403 (2002). 
[13] Y. Shin, M. Saba ,M. Vengalattore, T. A. Pasquini , C. 

Sanner, A. E. Leanhardt , M. Prentiss, D. E. Pritchard, 

and W. Ketterle, Phys. Rev. Lett. (2004) 160406. 
[14] A. Ferrando, M. Zacares, and M. A. Garcia-March, J. 

A. Monsoriu, and P. Fernandez de Cordoba, Phys. Rev. 

Lett. 95, 123901 (2005). 
[15] Y. V. Kartashov, A. Ferrando, A. A. Egorov, and LI. 

Torner, Phys. Rev. Lett. 95, 123902 (2005). 
[16] M. Oster, M. Johansson, Phys. Rev. E 73, 066608 (2006). 
[17] V. M. Perez-Garcia, H. Michinel, H. Herrero, Phys. Rev. 

A 57, 3837 (1998). 
[18] J. J. Garcia-Ripoll, G. Molina- Terriza, V. M. Perez- 

Garcia, and LI. Torner, Phys. Rev. Lett. 87, 140403 

(2001). 

[19] G. D. Montesinos, V. M. Perez-Garcia, Math. Comp. 
Simul. 69, 447-456 (2005). 



